ADNP dysregulates methylation and mitochondrial gene expression in the cerebellum of a Helsmoortel-Van der Aa syndrome autopsy case.

BACKGROUND: Helsmoortel-Van der Aa syndrome is a neurodevelopmental disorder in which patients present with autism, intellectual disability, and frequent extra-neurological features such as feeding and gastrointestinal problems, visual impairments, and cardiac abnormalities. All patients exhibit heterozygous de novo nonsense or frameshift stop mutations in the Activity-Dependent Neuroprotective Protein (ADNP) gene, accounting for a prevalence of 0.2% of all autism cases worldwide. ADNP fulfills an essential chromatin remodeling function during brain development. In this study, we investigated the cerebellum of a died 6-year-old male patient with the c.1676dupA/p.His559Glnfs*3 ADNP mutation. RESULTS: The clinical presentation of the patient was representative of the Helsmoortel-Van der Aa syndrome. During his lifespan, he underwent two liver transplantations after which the child died because of multiple organ failure. An autopsy was performed, and various tissue samples were taken for further analysis. We performed a molecular characterization of the cerebellum, a brain region involved in motor coordination, known for its highest ADNP expression and compared it to an age-matched control subject. Importantly, epigenome-wide analysis of the ADNP cerebellum identified CpG methylation differences and expression of multiple pathways causing neurodevelopmental delay. Interestingly, transcription factor motif enrichment analysis of differentially methylated genes showed that the ADNP binding motif was the most significantly enriched. RNA sequencing of the autopsy brain further identified downregulation of the WNT signaling pathway and autophagy defects as possible causes of neurodevelopmental delay. Ultimately, label-free quantification mass spectrometry identified differentially expressed proteins involved in mitochondrial stress and sirtuin signaling pathways amongst others. Protein-protein interaction analysis further revealed a network including chromatin remodelers (ADNP, SMARCC2, HDAC2 and YY1), autophagy-related proteins (LAMP1, BECN1 and LC3) as well as a key histone deacetylating enzyme SIRT1, involved in mitochondrial energy metabolism. The protein interaction of ADNP with SIRT1 was further biochemically validated through the microtubule-end binding proteins EB1/EB3 by direct co-immunoprecipitation in mouse cerebellum, suggesting important mito-epigenetic crosstalk between chromatin remodeling and mitochondrial energy metabolism linked to autophagy stress responses. This is further supported by mitochondrial activity assays and stainings in patient-derived fibroblasts which suggest mitochondrial dysfunctions in the ADNP deficient human brain. CONCLUSION: This study forms the baseline clinical and molecular characterization of an ADNP autopsy cerebellum, providing novel insights in the disease mechanisms of the Helsmoortel-Van der Aa syndrome. By combining multi-omic and biochemical approaches, we identified a novel SIRT1-EB1/EB3-ADNP protein complex which may contribute to autophagic flux alterations and impaired mitochondrial metabolism in the Helsmoortel-Van der Aa syndrome and holds promise as a new therapeutic target.

Loss-of-function of activity-dependent neuroprotective protein (ADNP) by a splice-acceptor site mutation causes Helsmoortel-Van der Aa syndrome.

Mutations in ADNP result in Helsmoortel-Van der Aa syndrome. Here, we describe the first de novo intronic deletion, affecting the splice-acceptor site of the first coding ADNP exon in a five-year-old girl with developmental delay and autism. Whereas exome sequencing failed to detect the non-coding deletion, genome-wide CpG methylation analysis revealed an episignature suggestive of a Helsmoortel-Van der Aa syndrome diagnosis. This diagnosis was further supported by PhenoScore, a novel facial recognition software package. Subsequent whole-genome sequencing resolved the three-base pair ADNP deletion c.[-5-1_-4del] with transcriptome sequencing showing this deletion leads to skipping of exon 4. An N-terminal truncated protein could not be detected in transfection experiments with a mutant expression vector in HEK293T cells, strongly suggesting this is a first confirmed diagnosis exclusively due to haploinsufficiency of the ADNP gene. Pathway analysis of the methylome indicated differentially methylated genes involved in brain development, the cytoskeleton, locomotion, behavior, and muscle development. Along the same line, transcriptome analysis identified most of the differentially expressed genes as upregulated, in line with the hypomethylated CpG episignature and confirmed the involvement of the cytoskeleton and muscle development pathways that are also affected in patient cell lines and animal models. In conclusion, this novel mutation for the first time demonstrates that Helsmoortel-Van der Aa syndrome can be caused by a loss-of-function mutation. Moreover, our study elegantly illustrates the use of EpiSignatures, WGS and Phenoscore as novel complementary diagnostic tools in case a of negative WES result.

Identification of a DLG3 stop mutation in the MRX20 family.

Here, we identified the causal mutation in the MRX20 family, one of the larger X-linked pedigrees that have been described in which no gene had been identified up till now. In 1995, the putative disease gene had been mapped to the pericentromeric region on the X chromosome, but no follow-up studies were performed. Here, whole exome sequencing (WES) on two affected and one unaffected family member revealed the c.195del/p.(Thr66ProfsTer55) mutation in the DLG3 gene (NM_021120.4) that segregated with the affected individuals in the family. DLG3 mutations have been consequently associated with intellectual disability and are a plausible explanation for the clinical abnormalities observed in this family. In addition, we identified two other variants co-segregating with the phenotype: a stop gain mutation in SSX1 (c.358G>T/p.(Glu120Ter)) (NM_001278691.2) and a nonsynonymous SNV in USP27X (c.56 A>G/p.(Gln19Arg)) (NM_001145073.3). RNA sequencing revealed 14 differentially expressed genes (p value < 0.1) in 7 affected males compared to 4 unaffected males of the family, including four genes known to be associated with neurological disorders. Thus, in this paper we identified the c.195del/p.(Thr66ProfsTer55) mutation in the DLG3 gene (NM_021120.4) as likely responsible for the phenotype observed in the MRX20 family.

An Additional Lrp4 High Bone Mass Mutation Mitigates the Sost-Knockout Phenotype in Mice by Increasing Bone Remodeling.

Pathogenic variants disrupting the binding between sclerostin (encoded by SOST) and its receptor LRP4 have previously been described to cause sclerosteosis, a rare high bone mass disorder. The sclerostin-LRP4 complex inhibits canonical WNT signaling, a key pathway regulating osteoblastic bone formation and a promising therapeutic target for common bone disorders, such as osteoporosis. In the current study, we crossed mice deficient for Sost (Sost-/-) with our p.Arg1170Gln Lrp4 knock-in (Lrp4KI/KI) mouse model to create double mutant Sost-/-;Lrp4KI/KI mice. We compared the phenotype of Sost-/- mice with that of Sost-/-;Lrp4KI/KI mice, to investigate a possible synergistic effect of the disease-causing p.Arg1170Trp variant in Lrp4 on Sost deficiency. Interestingly, presence of Lrp4KI alleles partially mitigated the Sost-/- phenotype. Cellular and dynamic histomorphometry did not reveal mechanistic insights into the observed phenotypic differences. We therefore determined the molecular effect of the Lrp4KI allele by performing bulk RNA sequencing on Lrp4KI/KI primary osteoblasts. Unexpectedly, mostly genes related to bone resorption or remodeling (Acp5, Rankl, Mmp9) were upregulated in Lrp4KI/KI primary osteoblasts. Verification of these markers in Lrp4KI/KI, Sost-/- and Sost-/-;Lrp4KI/KI mice revealed that sclerostin deficiency counteracts this Lrp4KI/KI effect in Sost-/-;Lrp4KI/KI mice. We therefore hypothesize that models with two inactivating Lrp4KI alleles rather activate bone remodeling, with a net gain in bone mass, whereas sclerostin deficiency has more robust anabolic effects on bone formation. Moreover, these effects of sclerostin and Lrp4 are stronger in female mice, contributing to a more severe phenotype than in males and more detectable phenotypic differences among different genotypes.

Mitochondrial GpC and CpG DNA Hypermethylation Cause Metabolic Stress-Induced Mitophagy and Cholestophagy.

Metabolic dysfunction-associated steatotic liver disease (MASLD) is characterized by a constant accumulation of lipids in the liver. This hepatic lipotoxicity is associated with a dysregulation of the first step in lipid catabolism, known as beta oxidation, which occurs in the mitochondrial matrix. Eventually, this dysregulation will lead to mitochondrial dysfunction. To evaluate the possible involvement of mitochondrial DNA methylation in this lipid metabolic dysfunction, we investigated the functional metabolic effects of mitochondrial overexpression of CpG (MSssI) and GpC (MCviPI) DNA methyltransferases in relation to gene expression and (mito)epigenetic signatures. Overall, the results show that mitochondrial GpC and, to a lesser extent, CpG methylation increase bile acid metabolic gene expression, inducing the onset of cholestasis through mito-nuclear epigenetic reprogramming. Moreover, both increase the expression of metabolic nuclear receptors and thereby induce basal overactivation of mitochondrial respiration. The latter promotes mitochondrial swelling, favoring lipid accumulation and metabolic-stress-induced mitophagy and autophagy stress responses. In conclusion, both mitochondrial GpC and CpG methylation create a metabolically challenging environment that induces mitochondrial dysfunction, which may contribute to the progression of MASLD.

Heterozygous variants in CTR9, which encodes a major component of the PAF1 complex, are associated with a neurodevelopmental disorder.

PURPOSE: CTR9 is a subunit of the PAF1 complex (PAF1C) that plays a crucial role in transcription regulation by binding CTR9 to RNA polymerase II. It is involved in transcription-coupled histone modification through promoting H3K4 and H3K36 methylation. We describe the clinical and molecular studies in 13 probands, harboring likely pathogenic CTR9 missense variants, collected through GeneMatcher. METHODS: Exome sequencing was performed in all individuals. CTR9 variants were assessed through 3-dimensional modeling of the activated human transcription complex Pol II-DSIF-PAF-SPT6 and the PAF1/CTR9 complex. H3K4/H3K36 methylation analysis, mitophagy assessment based on tetramethylrhodamine ethyl ester perchlorate immunofluorescence, and RNA-sequencing in skin fibroblasts from 4 patients was performed. RESULTS: Common clinical findings were variable degrees of intellectual disability, hypotonia, joint hyperlaxity, speech delay, coordination problems, tremor, and autism spectrum disorder. Mild dysmorphism and cardiac anomalies were less frequent. For 11 CTR9 variants, de novo occurrence was shown. Three-dimensional modeling predicted a likely disruptive effect of the variants on local CTR9 structure and protein interaction. Additional studies in fibroblasts did not unveil the downstream functional consequences of the identified variants. CONCLUSION: We describe a neurodevelopmental disorder caused by (mainly) de novo variants in CTR9, likely affecting PAF1C function.

Shared brain transcriptomic signature in TDP-43 type A FTLD patients with or without GRN mutations.

Frontotemporal lobar degeneration with TDP-43 inclusions (FTLD-TDP) is a complex heterogeneous neurodegenerative disorder for which mechanisms are poorly understood. To explore transcriptional changes underlying FTLD-TDP, we performed RNA-sequencing on 66 genetically unexplained FTLD-TDP patients, 24 FTLD-TDP patients with GRN mutations and 24 control participants. Using principal component analysis, hierarchical clustering, differential expression and coexpression network analyses, we showed that GRN mutation carriers and FTLD-TDP-A patients without a known mutation shared a common transcriptional signature that is independent of GRN loss-of-function. After combining both groups, differential expression as compared to the control group and coexpression analyses revealed alteration of processes related to immune response, synaptic transmission, RNA metabolism, angiogenesis and vesicle-mediated transport. Deconvolution of the data highlighted strong cellular alterations that were similar in FTLD-TDP-A and GRN mutation carriers with NSF as a potentially important player in both groups. We propose several potentially druggable pathways such as the GABAergic, GDNF and sphingolipid pathways. Our findings underline new disease mechanisms and strongly suggest that affected pathways in GRN mutation carriers extend beyond GRN and contribute to genetically unexplained forms of FTLD-TDP-A.

Induced pluripotent stem cell-derived motor neurons of CMT type 2 patients reveal progressive mitochondrial dysfunction.

Axonal Charcot-Marie-Tooth neuropathies (CMT type 2) are caused by inherited mutations in various genes functioning in different pathways. The types of genes and multiplicity of mutations reflect the clinical and genetic heterogeneity in CMT2 disease, which complicates its diagnosis and has inhibited the development of therapies. Here, we used CMT2 patient-derived pluripotent stem cells (iPSCs) to identify common hallmarks of axonal degeneration shared by different CMT2 subtypes. We compared the cellular phenotypes of neurons differentiated from CMT2 patient iPSCs with those from healthy controls and a CRISPR/Cas9-corrected isogenic line. Our results demonstrated neurite network alterations along with extracellular electrophysiological abnormalities in the differentiated motor neurons. Progressive deficits in mitochondrial and lysosomal trafficking, as well as in mitochondrial morphology, were observed in all CMT2 patient lines. Differentiation of the same CMT2 iPSC lines into peripheral sensory neurons only gave rise to cellular phenotypes in subtypes with sensory involvement, supporting the notion that some gene mutations predominantly affect motor neurons. We revealed a common mitochondrial dysfunction in CMT2-derived motor neurons, supported by alterations in the expression pattern and oxidative phosphorylation, which could be recapitulated in the sciatic nerve tissue of a symptomatic mouse model. Inhibition of a dual leucine zipper kinase could partially ameliorate the mitochondrial disease phenotypes in CMT2 subtypes. Altogether, our data reveal shared cellular phenotypes across different CMT2 subtypes and suggests that targeting such common pathomechanisms could allow the development of a uniform treatment for CMT2.

The role of endothelial autocrine NRG1/ERBB4 signaling in cardiac remodeling.

Neuregulin-1 (NRG1) is a paracrine growth factor, secreted by cardiac endothelial cells (ECs) in conditions of cardiac overload/injury. The current concept is that the cardiac effects of NRG1 are mediated by activation of erythroblastic leukemia viral oncogene homolog (ERBB)4/ERBB2 receptors on cardiomyocytes. However, recent studies have shown that paracrine effects of NRG1 on fibroblasts and macrophages are equally important. Here, we hypothesize that NRG1 autocrine signaling plays a role in cardiac remodeling. We generated EC-specific Erbb4 knockout mice to eliminate endothelial autocrine ERBB4 signaling without affecting paracrine NRG1/ERBB4 signaling in the heart. We first observed no basal cardiac phenotype in these mice up to 32 wk. We next studied these mice following transverse aortic constriction (TAC), exposure to angiotensin II (ANG II), or myocardial infarction in terms of cardiac performance, myocardial hypertrophy, myocardial fibrosis, and capillary density. In general, no major differences between EC-specific Erbb4 knockout mice and control littermates were observed. However, 8 wk following TAC both myocardial hypertrophy and fibrosis were attenuated by EC-specific Erbb4 deletion, albeit these responses were normalized after 20 wk. Similarly, 4 wk after ANG II treatment, myocardial fibrosis was less pronounced compared with control littermates. These observations were supported by RNA-sequencing experiments on cultured endothelial cells showing that NRG1 controls the expression of various hypertrophic and fibrotic pathways. Overall, this study shows a role of endothelial autocrine NRG1/ERBB4 signaling in the modulation of hypertrophic and fibrotic responses during early cardiac remodeling. This study contributes to understanding the spatiotemporal heterogeneity of myocardial autocrine and paracrine responses following cardiac injury.NEW & NOTEWORTHY The role of NRG1/ERBB signaling in endothelial cells is not completely understood. Our study contributes to the understanding of spatiotemporal heterogeneity of myocardial autocrine and paracrine responses following cardiac injury and shows a role of endothelial autocrine NRG1/ERBB4 signaling in the modulation of hypertrophic and fibrotic responses during early cardiac remodeling.

Transcriptional dysregulation by a nucleus-localized aminoacyl-tRNA synthetase associated with Charcot-Marie-Tooth neuropathy.

Charcot-Marie-Tooth disease (CMT) is a length-dependent peripheral neuropathy. The aminoacyl-tRNA synthetases constitute the largest protein family implicated in CMT. Aminoacyl-tRNA synthetases are predominantly cytoplasmic, but are also present in the nucleus. Here we show that a nuclear function of tyrosyl-tRNA synthetase (TyrRS) is implicated in a Drosophila model of CMT. CMT-causing mutations in TyrRS induce unique conformational changes, which confer capacity for aberrant interactions with transcriptional regulators in the nucleus, leading to transcription factor E2F1 hyperactivation. Using neuronal tissues, we reveal a broad transcriptional regulation network associated with wild-type TyrRS expression, which is disturbed when a CMT-mutant is expressed. Pharmacological inhibition of TyrRS nuclear entry with embelin reduces, whereas genetic nuclear exclusion of mutant TyrRS prevents hallmark phenotypes of CMT in the Drosophila model. These data highlight that this translation factor may contribute to transcriptional regulation in neurons, and suggest a therapeutic strategy for CMT.

Cardiac endothelial cell transcriptome in neonatal, adult, and remodeling hearts.

Cardiac microvascular endothelial cells (CMVECs) are the most numerous cells in the myocardium and orchestrate cardiogenesis during development, regulate adult cardiac function, and modulate pathophysiology of heart failure. It has been shown that the transcriptome of CMVECs differs from other endothelial cell types, but transcriptomic changes in cardiac endothelial cells during cardiac maturation and cardiac remodeling have not been studied. CMVECs were isolated from rat hearts based on CD31 expression and were immediately processed for RNA sequencing. We compared gene expression levels from primary CMVECs of neonatal hearts, normal adult hearts, and infarcted hearts. Between neonatal and adult CMVECs, 6,838 genes were differentially expressed, indicating that CMVECs undergo a substantial transformation during postnatal cardiac growth. A large fraction of genes upregulated in neonatal CMVECs are part of mitosis pathways, whereas a large fraction of genes upregulated in adult CMVECs are part of cellular response, secretory, signaling, and cell adhesion pathways. Between CMVECs of normal adult hearts and infarcted hearts, 159 genes were differentially expressed. We found a limited degree of overlap (55 genes) between the differentially expressed genes in neonatal and infarcted-hearts. Of 46 significantly upregulated genes in the infarcted heart, 46% were also upregulated in neonatal hearts relative to sham. Of 113 significantly downregulated genes in the infarcted-hearts, 30% were also downregulated in neonatal hearts relative to sham. These data demonstrate that CMVECs undergo dramatic changes from neonatal to adult and more subtle changes between normal state and cardiac remodeling.

Mutant HSPB1 causes loss of translational repression by binding to PCBP1, an RNA binding protein with a possible role in neurodegenerative disease.

The small heat shock protein HSPB1 (Hsp27) is an ubiquitously expressed molecular chaperone able to regulate various cellular functions like actin dynamics, oxidative stress regulation and anti-apoptosis. So far disease causing mutations in HSPB1 have been associated with neurodegenerative diseases such as distal hereditary motor neuropathy, Charcot-Marie-Tooth disease and amyotrophic lateral sclerosis. Most mutations in HSPB1 target its highly conserved α-crystallin domain, while other mutations affect the C- or N-terminal regions or its promotor. Mutations inside the α-crystallin domain have been shown to enhance the chaperone activity of HSPB1 and increase the binding to client proteins. However, the HSPB1-P182L mutation, located outside and downstream of the α-crystallin domain, behaves differently. This specific HSPB1 mutation results in a severe neuropathy phenotype affecting exclusively the motor neurons of the peripheral nervous system. We identified that the HSPB1-P182L mutant protein has a specifically increased interaction with the RNA binding protein poly(C)binding protein 1 (PCBP1) and results in a reduction of its translational repressive activity. RNA immunoprecipitation followed by RNA sequencing on mouse brain lead to the identification of PCBP1 mRNA targets. These targets contain larger 3'- and 5'-UTRs than average and are enriched in an RNA motif consisting of the CTCCTCCTCCTCC consensus sequence. Interestingly, next to the clear presence of neuronal transcripts among the identified PCBP1 targets we identified known genes associated with hereditary peripheral neuropathies and hereditary spastic paraplegias. We therefore conclude that HSPB1 can mediate translational repression through interaction with an RNA binding protein further supporting its role in neurodegenerative disease.

Charcot-Marie-Tooth disease type 2G redefined by a novel mutation in LRSAM1.

OBJECTIVE: To identify the unknown genetic cause in a large pedigree previously classified with a distinct form of axonal Charcot-Marie-Tooth disease type 2G (CMT2G) and to explore its transcriptional consequences. METHODS: Clinical reevaluation of the pedigree was performed, followed by linkage analysis with the redefined disease statuses, and whole genome and exome sequencing. The impact of the mutation was investigated by immunoblotting and transcriptome sequencing. RESULTS: Thirteen affected individuals over 3 generations displayed mild and quiescent lower-limb axonal sensorimotor neuropathy. Magnetic resonance imaging (MRI) of lower-limb musculature systematically showed fatty atrophy in clinical and subclinical mutation carriers. We redefined the disease-linked region to chr9q31.3-q34.2 and subsequently identified a novel missense variant in the E3 ubiquitin-protein ligase LRSAM1 (p.Cys694Tyr). Unlike previous reports, we demonstrated in patients' lymphoblasts that the mutation does not influence overall protein levels of LRSAM1, nor of its ubiquitylation target TSG101. The mutation is associated with several transcriptional changes, including a significant upregulation of another E3 ubiquitin-protein ligase, NEDD4L, and of TNFRSF21, a key regulator of axonal degeneration. INTERPRETATION: Our findings demonstrate that the isolated genetic entity CMT2G is caused by a missense mutation in LRSAM1 and should be reclassified as CMT2P. MRI of lower-limb musculature can be used to detect minimal signs of the disease. Transcriptome analysis of patients' cells highlights novel molecular players associated with LRSAM1 dysfunction, and reveals pathways and therapeutic targets shared with amyotrophic lateral sclerosis and Alzheimer disease. Ann Neurol 2016;80:823-833.

A comprehensive study of the genetic impact of rare variants in SORL1 in European early-onset Alzheimer's disease.

The sortilin-related receptor 1 (SORL1) gene has been associated with increased risk for Alzheimer's disease (AD). Rare genetic variants in the SORL1 gene have also been implicated in autosomal dominant early-onset AD (EOAD). Here we report a large-scale investigation of the contribution of genetic variability in SORL1 to EOAD in a European EOAD cohort. We performed massive parallel amplicon-based re-sequencing of the full coding region of SORL1 in 1255 EOAD patients and 1938 age- and origin-matched control individuals in the context of the European Early-Onset Dementia (EOD) consortium, originating from Belgium, Spain, Portugal, Italy, Sweden, Germany, and Czech Republic. We identified six frameshift variants and two nonsense variants that were exclusively present in patients. These mutations are predicted to result in haploinsufficiency through nonsense-mediated mRNA decay, which could be confirmed experimentally for SORL1 p.Gly447Argfs*22 observed in a Belgian EOAD patient. We observed a 1.5-fold enrichment of rare non-synonymous variants in patients (carrier frequency 8.8 %; SkatOMeta p value 0.0001). Of the 84 non-synonymous rare variants detected in the full patient/control cohort, 36 were only detected in patients. Our findings underscore a role of rare SORL1 variants in EOAD, but also show a non-negligible frequency of these variants in healthy individuals, necessitating the need for pathogenicity assays. Premature stop codons due to frameshift and nonsense variants, have so far exclusively been found in patients, and their predicted mode of action corresponds with evidence from in vitro functional studies of SORL1 in AD.

Concurrent whole-genome haplotyping and copy-number profiling of single cells.

Methods for haplotyping and DNA copy-number typing of single cells are paramount for studying genomic heterogeneity and enabling genetic diagnosis. Before analyzing the DNA of a single cell by microarray or next-generation sequencing, a whole-genome amplification (WGA) process is required, but it substantially distorts the frequency and composition of the cell's alleles. As a consequence, haplotyping methods suffer from error-prone discrete SNP genotypes (AA, AB, BB) and DNA copy-number profiling remains difficult because true DNA copy-number aberrations have to be discriminated from WGA artifacts. Here, we developed a single-cell genome analysis method that reconstructs genome-wide haplotype architectures as well as the copy-number and segregational origin of those haplotypes by employing phased parental genotypes and deciphering WGA-distorted SNP B-allele fractions via a process we coin haplarithmisis. We demonstrate that the method can be applied as a generic method for preimplantation genetic diagnosis on single cells biopsied from human embryos, enabling diagnosis of disease alleles genome wide as well as numerical and structural chromosomal anomalies. Moreover, meiotic segregation errors can be distinguished from mitotic ones.

Frequency and pattern of heteroplasmy in the complete human mitochondrial genome.

Determining the levels of human mitochondrial heteroplasmy is of utmost importance in several fields. In spite of this, there are currently few published works that have focused on this issue. In order to increase the knowledge of mitochondrial DNA (mtDNA) heteroplasmy, the main goal of this work is to investigate the frequency and the mutational spectrum of heteroplasmy in the human mtDNA genome. To address this, a set of nine primer pairs designed to avoid co-amplification of nuclear DNA (nDNA) sequences of mitochondrial origin (NUMTs) was used to amplify the mitochondrial genome in 101 individuals. The analysed individuals represent a collection with a balanced representation of genders and mtDNA haplogroup distribution, similar to that of a Western European population. The results show that the frequency of heteroplasmic individuals exceeds 61%. The frequency of point heteroplasmy is 28.7%, with a widespread distribution across the entire mtDNA. In addition, an excess of transitions in heteroplasmy were detected, suggesting that genetic drift and/or selection may be acting to reduce its frequency at population level. In fact, heteroplasmy at highly stable positions might have a greater impact on the viability of mitochondria, suggesting that purifying selection must be operating to prevent their fixation within individuals. This study analyses the frequency of heteroplasmy in a healthy population, carrying out an evolutionary analysis of the detected changes and providing a new perspective with important consequences in medical, evolutionary and forensic fields.

Single-cell paired-end genome sequencing reveals structural variation per cell cycle.

The nature and pace of genome mutation is largely unknown. Because standard methods sequence DNA from populations of cells, the genetic composition of individual cells is lost, de novo mutations in cells are concealed within the bulk signal and per cell cycle mutation rates and mechanisms remain elusive. Although single-cell genome analyses could resolve these problems, such analyses are error-prone because of whole-genome amplification (WGA) artefacts and are limited in the types of DNA mutation that can be discerned. We developed methods for paired-end sequence analysis of single-cell WGA products that enable (i) detecting multiple classes of DNA mutation, (ii) distinguishing DNA copy number changes from allelic WGA-amplification artefacts by the discovery of matching aberrantly mapping read pairs among the surfeit of paired-end WGA and mapping artefacts and (iii) delineating the break points and architecture of structural variants. By applying the methods, we capture DNA copy number changes acquired over one cell cycle in breast cancer cells and in blastomeres derived from a human zygote after in vitro fertilization. Furthermore, we were able to discover and fine-map a heritable inter-chromosomal rearrangement t(1;16)(p36;p12) by sequencing a single blastomere. The methods will expedite applications in basic genome research and provide a stepping stone to novel approaches for clinical genetic diagnosis.

Genome-wide copy number profiling of single cells in S-phase reveals DNA-replication domains.

Single-cell genomics is revolutionizing basic genome research and clinical genetic diagnosis. However, none of the current research or clinical methods for single-cell analysis distinguishes between the analysis of a cell in G1-, S- or G2/M-phase of the cell cycle. Here, we demonstrate by means of array comparative genomic hybridization that charting the DNA copy number landscape of a cell in S-phase requires conceptually different approaches to that of a cell in G1- or G2/M-phase. Remarkably, despite single-cell whole-genome amplification artifacts, the log2 intensity ratios of single S-phase cells oscillate according to early and late replication domains, which in turn leads to the detection of significantly more DNA imbalances when compared with a cell in G1- or G2/M-phase. Although these DNA imbalances may, on the one hand, be falsely interpreted as genuine structural aberrations in the S-phase cell's copy number profile and hence lead to misdiagnosis, on the other hand, the ability to detect replication domains genome wide in one cell has important applications in DNA-replication research. Genome-wide cell-type-specific early and late replicating domains have been identified by analyses of DNA from populations of cells, but cell-to-cell differences in DNA replication may be important in genome stability, disease aetiology and various other cellular processes.

Nuclear insertions of mitochondrial origin: Database updating and usefulness in cancer studies.

Nuclear insertions of mitochondrial origin (NUMTs) can be useful tools in evolution and population studies. However, due to their similarity to mitochondrial DNA (mtDNA), NUMTs may also be a source of contamination in mtDNA studies. The main goal of this work is to present a database of NUMTs, based on the latest version of the human genome-GRCh37 draft. A total of 755 insertions were identified. There are 33 paralogous sequences with over 80% sequence similarity and of a greater length than 500bp. The non-identical positions between paralogous sequences are listed for the first time. As an application example, the described database is used to evaluate the impact of NUMT contamination in cancer studies. The evaluation reveals that 220 positions from 256 with zero hits in the current mtDNA phylogeny could in fact be traced to one or more nuclear insertions of mtDNA. This is due to they are located in non-identical positions between mtDNA and nuclear DNA (nDNA). After in silico primer validation of each revised cancer study, risk of co-amplification between mtDNA and nDNA was detected in some cases, whereas in others no risk of amplification was identified. This approach to cancer studies clearly proves the potential of our NUMT database as a valuable new tool to validate mtDNA mutations described in different contexts. Moreover, due to the amount of information provided for each nuclear insertion, this database should play an important role in designing evolutionary, phylogenetic and epidemiological studies.

Validated primer set that prevents nuclear DNA sequences of mitochondrial origin co-amplification: a revision based on the New Human Genome Reference Sequence (GRCh37).

A new human genome reference sequence--GRCh37--was recently generated and made available by the Genome Reference Consortium. Since the prior disposable human reference sequence--hg18--was previously used for the mitochondrial DNA primer BLAST validation, a revision of those previously published primer pairs is required. Thus, the aim of this Short Communication is to perform an in silico BLAST test of the published disposable nine primer pairs using the new human reference sequence and to report the pertinent modifications. The new analysis showed that one of the tested primer pairs requires a revision. Therefore, a new validated primer pair, which specifically amplifies the mitochondrial region located between positions 6520 and 9184, is presented.